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Single-vertex Feynman diagrams represent the dominant contribution to physical processes, but 
are frequently forbidden kinematically. This is changed when the particles involved propagate in a 
gravitational background and acquire an effective mass. Procedures are introduced that allow the 
calculation of lowest order diagrams, their corresponding transition probabilities, emission powers 
and spectra to all orders in the metric deviation, for particles of any spin propagating in gravitational 
fields described by any metric. Physical properties of the ’’space-time medium” are also discussed. 
It is shown in particular that a small dissipation term in the particle wave equations can trigger 
a strong back-reaction that introduces resonances in the radiative process and affects the resulting 
gravitational background. 


INTRODUCTION 


The search for particle processes of astrophysical significance to which gravitation makes a non-negligible contribu¬ 
tion is in general difficult, though potentially rewarding. The entire field of gravitational lensing is in fact based on 
the discovery of one such process in which light interacts with gravity [l[. It is hoped that a number of additional 
processes will point out new directions of investigation. The search is made more difficult, unfortunately, by the fact 
that the lowest order Feynman diagrams that represent the dominant contribution to a process, are frequently forbid¬ 
den by kinematics. Consider, for instance, the process in which an incoming massive particle of momentum and 
dispersion relation = m\ produces a photon of momentum = 0, while the outgoing particle has momentum 
= m?. Conservation of energy-momentum requires = p'^ + ifj_. In the rest frame of p we have 0 = p' -I- 
which gives i = —p', po = mi and again mi = \/P + m^ -|- i. Then (mi — leads to £ = (m^ — m^)/2mi 

which shows that for mi = m, the case considered, we get £ = 0 and the process becomes physically meaningless. 
There are processes, however, in which massive particles emitting a photon are not kinematically forbidden. This 
is certainly the case when gravitation alters the dispersion relations of at least one of the particles involved Q. We 
stress here that even the reduction of a Feynman diagram by a single vertex would result in a cross-section gain of a 
factor {GM/RY, where M and R are mass and radius of the gravitational source (units ?i = c = 1). 

External gravitational fields have long been known to play the role of a medium (3|,1^ in the propagation of particles, 
be these treated classically, or according to quantum mechanics. In the latter case scattering by a Newtonian potential 
has been the subject of several investigations but bremsstrahlung iOl , the emission of Cerenkov radiation 

[l^ and other processes [l^ have also been studied in connection with various gravitational sources. As stated above, 
external gravitational fields do alter the dispersion relations of a particle propagating in a gravitational background at 
least to the extent that the particle is no longer on shell. It would therefore appear sufficient to somehow solve a wave 
equation to obtain the desired result. This is done in quantum electrodynamics where the electro mag netic field is 
however represented by a four-vector of which only a single component is usually taken into account |14l | . The case of 
gravity, represented by a second rank tensor, is considerably more complicated. Moreover, the theoretical prediction 
by Mashhoon confirmed by other authors l^l^, of the existence of a coupling of gravity to spin, requires 

that the effect of a gravitational field be no longer limited to a single component of the metric. This has become even 
more pressing since the experimental observation of spin-rotation coupling for photons and neutrons 241 and of 
other important spin-induced effects at the macroscopic level [2H - E^ . 

External gravitational fields contribute to the solution of covariant wave equations through a Berry phase 28, 2^. 
This should be expected because in metric theories of gravitation [s^, general relativity in particular, the space 
parameter of Berry’s theory coincides with space-time. It has been shown that the wave equations for fermions and 
bosons can be solved exactly to first order in the metric deviation 7^1/ = 5^1/ — 77^^ for any metric Q„,,y and that 
the phases so calculated [iol give reliable results in interferometry, gyroscopy [i^ and optics 35, 3^, give the 


correct Einstein deflection, can be used in the study of neutrino helicity and flavour oscillations 3711 and of spin-gravity 
coupling in general Si], . They also reproduce a variety of known effects as discussed in [3^ l33| and [38| . 

The dispersion relations of a particle propagating in a gravitational background can be derived from the respective 
covariant wave equations. The gravitational phases mentioned above, change, in effect, a particle four-momentum by 
acting on the wave function of the field-free equations. This result applies equally well to fermions and bosons and can 
be extended to all orders in 7^^, . The calculation of even the most elementary Feynman diagrams does require an 
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FIG. 1: p' and i are the outgoing fermion and photon and P indicates the incoming fermion. 


appropriate treatment when gravitational fields are present. The procedures developed in Q fill in part this gap and 
apply to linearized gravitational, or inertial fields of any type up to intermediate intensities. The present paper focuses 
on the applications of the procedures outlined in Q rather than on the individual reactions, though the aim still is to 
find physical processes capable of leading to potentially observable results. Below, we give examples of processes that 
can be treated in ways that are similar. Because of the intrinsic resemblance, gravitational bremsstrahlung [ii-il 
will be added to this category in due time. 

Additional interesting possibilities do moreover arise when particles propagate in a gravitational background. In 
the conclusions, for instance, we briefly discuss properties of the ’’space-time medium” such as dispersion. 


THE PROCESS P ^ p' +'y 


Let us assume, for simplicity, that P in Fig.l is an incoming fermion and that the photon i and outgoing fermion 
p' are produced on-shell. The solution of the covariant Dirac equation for P, exact to 0{'y^v)-i is 


321 


^( 2 ;) = - m) e 

2m 


-i^T 


«'o(x) = f ^-0 , 


( 1 ) 


where -|-ir^(a:), is the covariant derivative, r^(a;) the spin connection and the matrices ^^{x) satisfy the 

relations { 7 ''(a;), 7 '^(a:)} = 2g^^'^{x). Both r^(x) and 7 '^(x) can be obtained from the usual constant Dirac matrices 
by using the tetrad fields and the relations 


7'^(x)=e^(x)7“ = (^^ + h^(x))7^ 


r^(x) 





( 2 ) 


where (t“^ = §[ 7 “, 7^]- A semicolon and a comma are also used as alternative ways to indicate covariant and partial 
derivatives respectively. We use units = c = 1, the signature of 77 ^ 1 / is —2 , = ‘hs + ‘I’G) 

^^{x)=V J dz^Tx{z) = J dz^[ya\,piz) - 'yfsx,a{z)]L°‘^{z)+ ^ J dz^j^xk^^ , (3) 

L°‘^{z) = (x“ — z°‘)k^ — (x^ — z^)k°‘ and 'I'o(a;) satisfies the usual, fiat space-time Dirac equation. It is convenient to 
re-write © in the form t['(x) = g{x)exp{—ipx)uo{p), where 

9{x) = ^ + hl{x)p& + <^gA^)) + H , (4) 


and 


1 

2 


dz^Af.x^p 


7/3A 




:7ai^P 
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(5) 
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We claim that the transition amplitude for the process of Fig.l can be calculated by introducing the generalized 
four-momentum 


P^i=P^. + h'^Pa + =p^, + Pf,, 


( 6 ) 


for the incoming fermion, as (|1]) itself suggests. The part that contains the gravitational field is indicated by 
In (ini), and <I>g are quantities that must be calculated, once the metric is known. They are related to the 

Fourier transforms of the corresponding expressions that appear in ([2]), ([3]) and ([5]). is not on-shell. In fact 

P'^P^ = ml = m^ + 2 , (7) 

where p^p^ = rrP because p^ is the momentum of the free fermion represented by 'I'o(2^) ia (HJ. The transition 
amplitude is then 

= -iZer]^^''uoip')ef,(^x)-fog{\^)uo{p ), ( 8 ) 


where q = p — p' — i, represents the polarization of the photon, and Ze is the charge of the fermion. It may 
be argued that a transition amplitude must be added to dH) at 0 ( 7 ^j/) because the contraction in dH) is in 

general accomplished by means ofo^*^ and g{x) contains a part that is independent of The transition amplitude 
M', estimated in 0 , is given by p 




= -iZe-f>^’'{\^)uo(p')ef,,^x)l0UoiP) ■ 


(9) 


However contains the part p^ of dl]), that comes from but not the new part that contains the gravitational 

contribution due to the propagation of the fermion in the field of the source. To this process is indistinguishable 

from bremsstrahlung and is dealt with, more properly, in that context. It will not be discussed further in this work. 

The calculation now requires that a metric be selected. 

Let us consider the particular instance of a fermion that is propagating with momentum p^ = p, impact parameter 
b > R and X 2 = 0, from 2:3 = —00 toward a gravitational source of mass M and radius R placed at the origin 
and described W the metric 700 = ‘^4> Tlij = where (j) = This metric is frequently used in leasing 


problems [44 1 


= (! + </.) 5'. All 


One finds Fq = — 1/2^,^ cr®-^ , F^ = — l/ 20 ,jtTb gO _ q ^ gO _ i_ 

spin matrices are now expressed in terms of ordinary, constant Dirac matrices. We also assume that the on-shell 
conditions p'^p'^ = m? = 0 remain valid. Extension of the calculation to include different particles, or higher 

order gravitational contributions to p' ,£ and (HD can be derived to all orders in [^. The Fourier transforms of 
the quantities that appear in (|4|) must now be calculated. We obtain 


KiPPa = 8Tr‘^6{qo)6{qx)6{qy)poGMKo{bqz), h^{q)pa = 8Tr‘^6{qo)S{qx)6{qy)pGMKoibq^). 


( 10 ) 


$G,2(g) = 0 , $G.3(g) = -8TT'^S{qo)S{q^)S{qy) -t p^ GMKoibq^). (11) 

Four-momentum conservation to zeroth order only is required because (ITO)) . (fTTI) are already of 0{^yv)- Though the 
gravitational field selected above is stationary, energy conservation must be introduced because the energy contribution 
of the field is contained in the generalized momentum of P. We further approximate the Bessel function Kglbqz) — 
x/Trj2bqze~^'^^ [1 — 1 / 86^2 -I- ...], itself a distribution, by Ko{bqz)/TT ~ 2Tr6{bqz) and eliminate S'^{q) from (fTUl) and (fTTI) . 
Conservation of energy-momentum will reappear as a factor (2TTpS^{q) in the expression for the radiated power W 
defined below. By removing (27r)"‘5^(g) from /i“ and 4 >g,^( 5') we obtain and $g.m of ([6]). We find 


GM 


Fb — po H- — Po + Pq ^ Pi — ^ 


GM (pI 


P 


p) =Pi,P2 = 0,P^ = P= p-(^p = p -I- P. 


GM 


( 12 ) 


We calculate the power radiated as photons in the process of Fig.l according to the formula [4£ 


There are two ways to calculate \Mf^p^\^. In the first one we replace Pa with P^ in the field-free = 0) expression 
given by E|Mp = 4 to^(p(jP“) -|- 8 (pap“)]. The gravitational contribution to M then appears in Py exclusively. 
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We also remove the terms —32m?(p'^p°‘) + 64m^ that do not contain gravitational contributions and therefore refer 
to the kinematically forbidden transition. This yields, to 0 ( 7 ^;/), the expression 




-4(p:,P“) + 8(p„P“) 


(14) 


In a second, alternate approach, we calculate |Mp directly ( 7 ^ 1 / 0) from ([ 8 |). By summing over final spins and 

averaging over initial spins and polarizations, we obtain 


= 2 H^'^^{(p' + ’^)7/3 (p + P + to) (^(p + m)^ + (p + m)P* + 77 (p + to)) 7 ^| , (15) 


2 ( 2 to ) 


where a = 7 ^^ and H = p^/p^(j){'y^p^ — j^p^) + (p°/p^ 7 ^p° — 'y^p^)4GMKi(bgz). On carrying out the traces of 
the Dirac matrices, the contribution from H vanishes. By further eliminating from |Mp the terms that refer to the 
kinematically forbidden transition, we again find (1141) . This supports our claim that the generalized momentum 
introduced in ([ 6 ]) leads to the correct value of the transition probability by the substitution of p^ with in the 
field-free expression. The integration over cPp' in (1131) is performed by means of the identity 




(16) 


while that over 0, the angle between P and £, can be carried out by writing the on-shell condition for p' in the form 

d{2\P\\£\cos9 - P°‘Pa + 2Poio + m^). (17) 

We find 


The radiation spectrum is given by 

Equation (113 and the condition 

I COS0| 


An \ b J p^ 


dW _ Z^£ (GM\ pl+p- 
d£ 2tt 


1 (2GM 


2pi^l - ^ I b 


{pl+p^)-2fpo (^1 + ^)|< 1 , 


(18) 


(19) 


( 20 ) 


lead, for p > to, to the inequality p{GM/b) < £ < p. It then follows that the hardest photons are emitted in the 
backward direction with energy £ ^ p and power 




2Z'^e^p^ 

Stt 



( 21 ) 


which would obviously take its highest values in the neighborhood of a very compact source. 
For p < m the inequality (1^ is satisfied for £{Gm/b) > to(1 — p/m) and we also find 




fGM\ 
47r \ 5 y 


pz y j 


( 22 ) 


which diverges for small values of p. This infrared divergence is well-known and arises as a consequence of the finite 
energy resolution Ae of the outgoing fermion. The process, as calculated, is in fact indistinguishable from that in 
which massless particles with energy < Ae are also emitted and from processes in which vertex corrections are present 
(virtual massless particles emitted and reabsorbed by the external lines of Fig.l). When these additional diagrams 
are calculated all infrared divergences disappear . In the particular case at hand, p in (l22l) can be simply replaced 
by GMpo/b. Below this value the process is not kinematically allowed. 

The process discussed in this section may be considered as the decay of a fermion of effective mass PaP°‘ into a 
photon and a fermion of mass to with a lifetime r = m/W which is indeed small. Quantitatively, for electrons with 
p ^ IGeV, £ ^ p, in the neighborhood of a canonical neutron star, we find r ~ 3 • 10“^®s. 
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THE PROCESS 7 ^ / + / 


Using the replacements I —>■ —L and P —>■ —q in Fig.l, we can calculate the process by which a photon produces 
a fermion / of momentum p' and an anti-fermion / of momentum q after propagating in a gravitational field. In 
addition —>■ —Z^e^ because of the presence of the antiparticle. Conservation of energy-momentum now requires 

that Lfj, = p'fi + q^: while the dispersion relations are p'^p'^ = qaq°‘ = and the generalized photon momentum is 
La = ia+^G,(T- In the center of mass frame of the (/, /)-system we now havep'-bg = L — Q and also Li = = 0. 

This and = 0 again show that if the effect of the gravitational field vanished = 0)) we would re-obtain the 
meaningless result £0 = -^ = 0. It is therefore the presence of the gravitational field that enables the process. The effect 
of the gravitational field on the polarization of the photon is given by Ea = -\- ^{^aa,i 3 — Li 3 < 7 ,a)S°'^ + § 7 a/ 3 ,cr'r“^ 

where the matrices 5'“^ and T°‘^ are given in and act on the matrix Oq,. The contributions of these terms reduce 
to \{d^a 2 — aJai)( 7 io -,2 — 72 cr.i)- In the approximation Ko{bK)/TT —>■ {2Tr)S{bK) to be used below, the derivatives of 
behave as {bK)6{bK) = 0 and the whole contribution of 5”“^ and to the photon polarization can be neglected. 
We assume that = 1^2 = 0, £3 = £ and use the following result 


$G.1 = l r d^°700.i^ r = 2GmUi + -j=^=) 

^ 7-00 ^ 7-00 0 V 0^ + 


which, for r ^ b, becomes $g,i 2GM£/b. Similarly, we find 

dz 


$G 3 = „ - GM£ 


^/WT. 


2GMi 


(&2 + z2)3/2 


2GM 


(23) 


(24) 


in the same approximation. Factorizing (27r)(5(&K), we find the generalized momenta Lq = £,Li ^ 2GM£/b, L 2 = 
0,^3= 7(1 — 2GM/b). By applying the substitutions indicated above, we can derive the transition amplitude for the 
process 7 —>■ / -b / 


= Z^e^ [4{p'^q^) + 8(g„g“)] , (25) 

and the rate at which energy is radiated as an anti-fermion 

"'■ = 8(S7 / SS'*'’ + ^ (26) 

The integration over d^p' can be carried out by means of the identity (1161) and the on-shell condition for p' becomes 

6 {{La - qa) (U“ - 9 “) - rn?) = S{LaL° - 2 Lo 9 ° + 2|T||9l cos^i). (27) 


Integrating over Oi, the angle between q and L, we find 

VFi = ^ / dq^-^ (l-LaL‘^ + . ( 28 ) 

47r J ^Lo\L\ \2 J 

We also hnd LaL°‘ ~ 27“$G.a —2GM£/b after replacing K[){bK)/'K with (27r)(5(6K), factorizing {2TT)5{bK) and writing 

1/|L| ~ (1 -b GMIb)l£. We finally obtain 


IFi 


87r 


GM\ m^q^ / 7^ 

~ J V ^ 


(29) 


from which can be immediately obtained. We also have 0(To — 9o) = 7 — 90 > 0, while | cos6>i| < 1 leads to 
0 < _(_ 2qQ which is always satisfied in the interval 7r/2 < < tt. 


THE PROCESS / + / ^ 7 

We now consider the process in which a fermion-antifermion couple in the initial state annihilates into a photon. 
We assume that the fermion that propagates in the gravitational background has generalized momentum P. By 
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conservation of energy-momentum we then have P^, = ip, and the generalized momentum is given by Pq = 

Po + poGM/b = po + Pq ,Pi = P 2 = 0 , P 3 = ^ 3(1 -I- GM/h) — p^GM/b = p 3 + P 3 . The transition amplitude becomes 


= ZV 4 ((,<,P") + 8 (p„P“ 


from which all terms referring to the kinematically forbidden transition have been eliminated. We also find 

o2 


W 2 = 


Z^e- 


d?qdH_^^^p ^ ^ f4(g„P“) + 8(p„P“) 


(27r)2 J qgPo 

The integration over d^q can be easily performed by using the identity J = J d^q&(qo)S(q^ — m^). We find 

W2 = J dHOiPo - io)^S {{Pa, - - n - m^) [po(r - P“) + 2(p„P“)] . 


(30) 


(31) 


The on-shell condition for q becomes 


5{q^ -m‘^) = 5 {{Pa, - ia,){P°‘ - n - m^) = I cos 62 + 


I J „ PaP“ - 24P0 - 


2l\P\ 


and the integration over 82 then yields 


W 2 = 


GM 
Att b 


di- 


P^\P\ 


{psPo + . 


In order to carry out the integration over £, we first calculate 


1 


1 


1 - 


Po|P| PsPo 

which must be substituted in (1341) . The integration over £ gives 

/GM 


1 /2GM Po GM 
2 \b P 3 b 


W 2 = 


Stt 


from which we obtain the radiation spectrum 

dW 2 Z^e^ [GM 


d£ 47r 

The condition | cos 02 1 < 1 then requires that 

{m^+P3Po)^ 


1 — m 


1 -m^ 


1 1 


P3P0 P 3 


1 1 


P3PQ PsJi 




Po + P3 + 


P3 b 


P 0 -P 3 + 


P3 b 


(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


Notice that the process calculated in this section is not the time-reversed of 7 —^ / -I- / because in the latter process 
gravitation is assumed to act on the incoming photon line and not on the outgoing fermion line. 


CONCLUSIONS 


External gravitational fields in radiative processes can be included in the calculation of a transition probability by 
simply replacing the momentum pp of a particle with its generalized version Pp in the corresponding expression for the 
zero-field process. The examples given involve spin-i and spin-1 particles, but the procedure can be extended to any 
spin. An essential point here is that the dispersion relations are altered by the extern al grav itational field and can be 
calculated if the corresponding wave equations can be solved to 0 ( 7 ^ 44 ), or higher [iSI l34l437j . It follows, in particular. 


that kinematically forbidden processes similar to that of FigH] become physical and their transition probabilities can 
be determined. The calculation of the gravitational contributions are greatly simplified and can be extended to higher 
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order in 7 ^ 1 ^. The applications are not confined to fields of a Newtonian type, but extend to any gravitational field. In 
this respect, the procedure presented goes beyond the results that apply to external electromagnetic potentials [l4l| . 


not only because the metric has in general ten components rather than just four, but also because time-independence 
is not a requirement. 

The transition amplitudes derived are to leading order and can therefore be considerably larger than those 

normally studied in the literature. The examples given are limited to a single type of physically relevant metric and 
we cannot conclude that the resulting spectra are general and can be used to identify the processes. The actual 
determination of the spectra requires the use of metrics specific to the problems studied. However, the results 
suggest that particle processes like bremsstrahlung, Cerenkov radiation, or positron production in the neighborhood 
of compact astrophysical objects, or in cosmology, need to be reconsidered. 

Space-time has so far been treated as a linear optical medium, though it is by no means clear what its ultimate 
properties will be as a result of quantum gravity. It is not in particular known whether its index of refraction will 
remain unaltered in response to high intensity fluxes of particles. There is scope for research on this and other 
properties of space-time. 

Our final considerations regard the back-reaction that physical processes may have on the gravitational background. 
We show below that the back reaction is not always negligible and provide an example of how a very small disturbance 
in the wave equation can grow rapidly and alter the background gravitational field. 

Equation o requires that 4*0 (a^) be a solution of the field-free Dirac equation and, of course, of the equation 


equations 


I^o(a:) = 0. The approximation procedure still holds true, however, when 'I'o(a;) satisfies more general 


With the addition of a dissipation term, the equation for 'Fq becomes 

+ ni^ - 2mado) 'Fq = 0 , 


(39) 


where we take a = a|('ko|7’|'I'o)r = [2l| and a is a dimensionless, arbitrary parameter, 0 < a < 1, that 


reflects the coupling strength of the dissipation term. When we substitute 'I'o(a:) = e^'pimaxo)(po{x) into (IMl) . we 
obtain 


[dl-dl+m^{l-a'^)\(l)o{x)=Q. (40) 

An example of a problem with similar behaviour is offered by a fluid heated from below. For small temperature 

gradients the fluid conducts the heat, but as the gradient increases conduction is not sufficient to lead the heat away 

and the fluid starts to convect. In realistic problems the exponential growth of 'I'q does not continue indefinitely, but 

is restricted at times Xq > t = lima by nonlinearities or dispersive effects that may have been initially neglected. 

The effect of the new solution vko on W can be found as follows. We first neglect the change m —>■ in W 

because in general cr < 1. Setting m = 1 for simplicity, the effect of the exponentially increasing term on m then 

2 

amounts to the transformations 5{Po -Po - 4) S{-2ia + Po-Po- io) and ^ (pokier - i 32 [(p‘i- 4 a^) 2 + 4 p‘ia '^] ) 
where we have used the relation /3 = p/po- W has therefore a resonance at po = 2a of width dtr^. Over times xq > t = 
{m/po)~^{GM/2b)~'^GeV~^, 4/ = T'^/q increases exponentially until the compensating mechanisms mentioned 
above kick in. For a proton of energy pq ^ lOGeF in the field of a canonical neutron star r ^ 3.5 x s. 

Considerably higher values of r can, of course, be obtained for the lighter fermions. As 4* grows, so does the energy 
momentum tensor associated with it and the gravitational field it generates, altering, in the process, the gravitational 
background. 
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